Supercurrent through a multi-level quantum dot close to singlet-triplet degeneracy 
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We investigate two serially-aligned quantum dots in the molecular regime of large tunnel couplings 
t. A Zeeman field B is used to tune the energy difference of singlet and triplet spin configurations. 
Attaching this geometry to BCS source and drain leads with gap A and phase difference gives rise 
to an equilibrium supercurrent J. To compute J in presence of Coulomb interactions U between 
the dot electrons, we employ the functional renormalization group (FRG). For B ~ t - where the 
singlet and (one out of a) triplet spin states are equal in energy - the current exhibits characteristics 
of a O-TT transition similar to a single impurity. Its magnitude in the n phase, however, jumps 
discontinuously at B = t, being smaller on the triplet side. Exploiting the flexibility of the FRG, we 
demonstrate that this effect is generic and calculate J for realistic experimental parameters A, U, 
and gate voltages e. To obtain a more thorough understanding of the discontinuity, we analytically 
treat the limit A = oo where one can access the exact many-particle states. Finally, carrying out 
perturbation theory in the dot-lead couplings substantiates the intuitive picture that Cooper pair 
tunneling is favored by a singlet spin configuration while inhibited by a triplet one. 

PACS numbers: 74.50. +r, 75.20.Hr 



I. INTRODUCTION 

Many intriguing phenomena in condensed matter 
physics are to be credited to many-body interactions. 
Two prominent examples are the Kondo effect as well 
as superconductivity. If a magnetic impurity is placed 
inside a superconducting metal, the low-energy behavior 
is governed by a competition of both. Depending on the 
ratio of the Kondo temperature Tk and the BCS gap 
A, Cooper pairs are broken up in order to screen the lo- 
cal spin {Tk ^ A, '0 phase'); or bulk superconductivity 
prevails and the impurity harbors a free magnetic mo- 
ment {Tk ^ A, 'tt phase'). At zero temperature T ~ 
and Zeeman field B = 0, a first order quantum phase 
transition from a non-degenerate to a twofold degenerate 
ground state takes place as the ratio A/Tk increases. 
This overall picture was essentially established decades 
ago^"— yet, the advent of nanotechnology now renders 
it feasible to fabricate quantum dots ('magnetic impu- 
rities') attached to BCS source and drain leads. The 
various microscopic parameters of these systems in ad- 
dition to Tk and A, such as a gate voltage e or left- 
right asymmetry in the dot-lead couplings F^ a, can be 
accurately determined (and possibly varied) in the ex- 
periments. This allows to study the physics - and thus 
correlations - in a very controlled way4i^— 

Consequently, a renewed interest in the interplay of 
superconductivity and the Kondo effect and how it man- 
ifests in the equilibrium supercurrent J (that generally 
flows between source and drain in presence of a finite 
phase difference (f>) for the particular quantum dot setup 
was triggered. Most theoretical studies focused on a sin- 
gle Anderson impurity coupled to BCS leads as a mini- 
mal (yet challenging) model,—"— and eventually both the 
boundary of the O-tt phase transition as well as J were cal- 
culated using reliable many-particle methods i^ "^"'^^" — 



By adjusting all microscopic parameters to the ones as- 
sociated with an experimental device, the critical super- 
current and hallmarks of the O-tt transition were indeed 
calculated in qualitative agreement with the measured 
data^i^ The functional renormalization group (FRG), 
which allows to flexibly treat arbitrary gate voltages e, 
gaps A, or asymmetries in the source-drain couplings 
^L,R, provided a convenient tool to achieve this goal. 

A realistic experimental quantum dot setup certainly 
features more than one single-particle levelj^ However, 
if the energy spacing S between subsequent levels is large 
compared to their broadening F, it is reasonable to as- 
sume that only the one closest to the Fermi energy de- 
termines the low-energy behavior. The experiments men- 
tioned above show a series of well-separated resonances 
and can apparently be described by such a single-level 
scenario; but the progress in measurement and device 
fabrication technologies makes it likely that multi-level 
quantum dots (in the sense that (5 < F) or more complex 
geometries where several dots are coupled will be subject 
of upcoming studies. This motivates us to address the 
intriguing physics of such multi-level setups with BCS 
leads theoretically. To this end, one needs to resort to 
multi-impurity Anderson modelS)^"— which over the last 
years were widely investigated for normal leads. 

In this paper, we study a geometry of two single-level 
dots aligned in series and attached to BCS leads. We 
focus on the parameter regime where the hopping t be- 
tween the dots is strong compared to the lead coupling 
^ F and explicitly demonstrate that the physics - e.g., a 
O-TT phase transition governing the Josephson current ~ 
can indeed be explained by single-level scenarios associ- 
ated with either the bonding or anti-bonding molecular 
states. One might naively hope that the same simple pic- 
ture emerges in presence of a Zeeman field B k, t where 
the single-particle (interaction C/ = 0) energies of bond- 
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FIG. 1. The multi-level quantum dot geometry subject to this 
paper. 



ing spin-up and anti-bonding spin-down states are equal; 
or rephrased, where the smallest two-particle eigenener- 
gies of an isolated dot (with U > 0), which are a spin 
singlet as well as one out of a triplet, are almost de- 
generate. However, this turns out to be only partially 
true: Whereas for i3 ss t and half filling various char- 
acteristics - e.g., the very idea of a level crossing phase 
transition as well as the corresponding lincshapes and 
parameter dependencies of the current - are just as they 
are in a single-level case, the magnitude of J changes 
discontinuously at -B = i in one of the phases. This indi- 
cates an additional first order (singlet-triplet) quantum 
phase transition. Treating the problem within a func- 
tional renormalization group approach, we illustrate that 
this scenario is generic and compute J as a function of 
the gate voltage for experimentally relevant values of U, 
A, and F. To gain a more thorough understanding of the 
singlet-triplet transition (which we expect to be a dis- 
tinct feature of any more complex quantum dot geome- 
try close a singlet-triplet degeneracy), we supplement our 
FRG analysis by an exact solution of the limit A = oo 
as well as with perturbation theory in the dot-lead cou- 
plings. 

The serial quantum dot geometry with BCS leads was 
recently studied in Rcf. '3l' using the numerical renormal- 
ization group and considering the limit t <^ T opposite 
to our one. Furthermore, Ref. [3l| only focused on the 
limit U = oo, B = 0, and fixed e. The idea of a singlet- 
triplet transition was discussed in Ref. |3^ for a general 
two-level geometry. However, a detailed investigation of 
this phenomenon, of how it manifests in J(e), and of how 
it can be interpreted intuitively in the context of a serial 
quantum dot for experimentally motivated values of U, 
A, and F is still missing. This work aims at filling this 
gap. 

We organize our paper as follows. The model is in- 
troduced in Sec. Ill Al The details of the functional RG 
approach, the treatment of the limit A = oo as well as 
perturbation theory in F are presented in Secs. lHI3llHC| 
and lHDl respectively. We present and discuss our results 
in Sec. IIHI and conclude with a brief summary. Finite- 
temperature data obtained at A = oo is shown in the 
Appendix. 



II. MODEL & METHOD 

A. Model 

Our multi-level quantum dot Josephson geometry is 
depicted in Fig. [TJ The BCS source and drain (i.e., left 



and right) leads are modeled by the Hamiltonian (in mo- 
mentum space and using standard second-quantized no- 
tation) 

k(T k 

(1) 
with A and (t)L,R = ^4>l'^ being the BCS gap and phase, 
respectively. The two Anderson impurities coupled by a 
hopping element t are described by 

i^dot = 5Z [(^ + ^) '^'t + (e - ^) n^x\ 

i=l,2 



(2) 



where nia = dj^d^^, and i = 1,2 labels a Wannier ba- 
sis. The gate voltage e = e — 3C//2 is shifted such that 
6 = corresponds to the point of particle-hole symmetry 
at zero Zeeman field B = 0.— Moreover, U denotes the 
strength of the Coulomb interaction. We refrain from 
introducing a level splitting and take local as well as 
nearest-neighbor interactions to be equal - yet not for 
fundamental reasons, but merely because it does not al- 
ter our main results. If necessary, both can be incorpo- 
rated with the functional RG approach. For the same 
reason, the dot is assumed to couple equally to the left 
and right lead (of size A^ — > oo; t is the local hopping 
strength): 
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A key ingredient to the methods discussed in the next 
Section is the dot Matsubara Green function QQiioj) as- 
sociated with the total Hamiltonian 



H = H, 
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in absence of Coulomb interactions {U ~ 0). It can be 
derived straightforwardly by expressing H in terms of 
Nambu spinors and subsequently employing equation of 
motion tcchniquesi^ The result reads 
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Assuming a constant local density of states p in the leads 
in absence of superconductivity (wide-band limit) gives 
rise to an energy-independent hybridization 
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B. Functional renormalization group 



A recently developed approach to compute low-energy 
transport properties of correlated quantum impurity sys- 
tems is provided by the functional renormalization group 
(FRG)." The latter reformulates a given many-particle 
problem in terms of an infinite set of coupled flow equa- 
tions for imaginary frequency (or time) single-particle ir- 
reducible vertex functions with an infrared cutoff as the 
flow parameter. Truncation of this hierarchy renders the 
FRG approximate w.r.t. the two-particle interaction and 
can hence a priori be justified only in the limit of small 
U . The most simple truncation scheme keeps track of the 
self-energy E'^ as well as of an effective Coulomb interac- 
tion U^ (i.e., the static part of the two-particle vertex). 
It yields flow equations for effective system parameters 
and can thus be regarded as an 'RG enhanced Hartree- 
Fock' approach - that does, however, not suffer from any 
mean field artifacts. Application of this scheme to various 
quantum dot geometries in equilibrium turned out to give 

I 



accurate results for the zero-temperature linear conduc- 
tance even at fairly large Coulomb interactions.— Finite 
temperatures or non-linear finite-bias transport can be 
addressed by employing a more elaborate truncations^— 



(7) or by a generahzation to Keldysh space^ 
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respec- 



tively. 

In the context of superconducting leads, the Hartree- 
Fock-like FRG approximation was shown to describe the 
T = physics - such as the boundary between the and 
■K phases or the Josephson current - of the single impu- 
rity case quantitatively at small to intermediate values 
of the Coulomb interaction and at least qualitatively at 
large U , again without being plagued by artifacts such 
as the breaking of spin symmetryi^ it is thus reasonable 
that the same holds for the present problem. We will 
a posteriori substantiate this by comparing with exact 
results obtained at A = oo. 

The zero-temperature FRG flow equations for the 
multi-level quantum dot geometry can be obtained by a 
straightforward generalization of the single-level case.— 
They read 
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for the self-energy as well as 
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for the effective Coulomb interaction. The subscripts j = 
1 ... 4 label the four Nambu indices associated with the 
electron spin and the two dot orbital degrees of freedom; 
the Green function at scale A is introduced as 



g^{iuj) = [go{iuj)-^-^ 



Ai-l 



(10) 



Supplemented by the initial conditions^i j]^^°° — 
as well as t/j2i2°° = ^2323°° = ~^1313°° = ^1414°° = 
^2323°° = ~^2424°° = —U (and all antisymmetric per- 
mutations), the ordinary coupled differential equations 
([8]) and ([9]) can be easily solved down to A = using 
standard numerical routines. 

The Josephson current is defined as J = — i([iJ, TV^J), 
with Nn being the particle number operator of the right 
lead; we take ft = 1 as well as the electron charge e = 1. 
It can be expressed analytically in terms of the exact 
interacting Green function Q{iuj)>^ 
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for which we eventually plug in the FRG approximation 



C. Exact solution at A = 00 

The limit of large BCS gaps allows for an exact 
solution.— i^i^i^ Namely, the non-interacting Green 
function of Eq. ([5]) becomes frequency-independent at 
A = 00 except for the trivial contribution itu, and one 
can therefore equivalently solve the Hamiltonian 
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dot 



e-^/^di^du 
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H.c. 



This can be achieved straightforwardly by again intro- 
ducing Nambu spinors, decomposing the 16-dimcnsional 
Hilbcrt space that underlies Hcb into sectors with dif- 
ferent Nambu particle numbers and subsequently com- 
puting all eigenvalues Ei analytically (the corresponding 
expressions are lengthy and will not be presented here). 
Moreover, it will prove helpful - particularly when sub- 
stantiating our interpretation of a singlet-triplet transi- 
tion - to characterize the associated eigenstates in terms 
of the square s ^ and the z-component s^ of the total spin: 
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where a = {ax, cTy, a^)^ is the usual vector containing the 
Pauh matrices o^^y^z- It is easy to show that s^ and s^ 
commute with each other as well as with /fefi; moreover, 
one should note that 2sz + 2 is just the Nambu particle 
number. Having calculated the set of eigenvalues Ei, the 
supercurrent can eventually be obtained from the phase 
derivative of the free energy>2i 
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D. Perturbation theory in F 

The functional renormalization group approximation 
is strictly controlled in the limit of small Coulomb inter- 
actions. A complementary approach (yet not RG-based) 
is to carry out an expansion w.r.t. the tunnel Hamilto- 
nian Ht, which can be justified if both the gap A as well 
as the distance between the many-particle eigcnstates of 
iJdot are large compared to rJ^i^i The resulting pertur- 
bation theory expression for the Josephson current can 
be interpreted in terms of successive tunnel events. Such 
intuitive picture will turn out to be helpful for under- 
standing the FRG results for the problem at hand. 

The current of Eg. pT|) can be easily expanded 
w.r.t. Ht (see Ref. [2l| for an instructive and detailed 
treatment of the single- level case). The first non- 
vanishing term reads 

J=-Im2j/ dTidT2dTz 
•3 ^ Jq 

X {THt{ti)Ht{t2)Ht{t^)H^^_)^ . 

(15) 

One can show that the imaginary part of the zcro- 
I 
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coupling thermal expectation value {. . . Jrl^'^")o is non- 
vanishing only if it additionally involves each of the op- 
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erators i^x.+j -f^T,+ 



H^'^_ once; subsequently plugging 



in the lead contribution {Tc\_^At 
out the wide-band limit yields 
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where E{e) = \J(^ + A^, and 
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We have simplified the last line assuming that T <C A 
but emphasize that one cannot neglect exp(— S/T) even 
in the zero-temperature limit as the imaginary time in- 
tegrations may give rise to a factor exp(£'/r). 

In order to evaluate the dot contribution 
<7dot('''i, ^2, Ts); which is complicated by the lack of 
Wick's theorem, one first solves the many-particle 
eigenvalue problem -ffdotl^) = A„|n). As we have 
assumed equal local- and nearest neighbor interac- 
tions, this can be achieved conveniently by rotating 
to a basis of bonding (b) and anti-bonding (a) states 
rfi,2CT = (rffccr ± daa)/^- After inserting unit operators 
1 = ^j^ l'T^)('^|i one can then perform the r-integrals for 
a given time ordering and eventually set T = 0. For 
Ti > T2 > T3 we obtain (and likewise for the other five 
possibilities) 



r0 rP rP i 
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[E{ei) + Xp\[E{ei) + A.][^(ei) + E{e2) + A„] [^(ea) + \p][E{e2) + \z\[E{ei) + ^(ea) + A,] 

I 

with A„ = A„ — Aq denoting the difference to the smallest intermediate states \n,p,q), a factor exp(A„_q_p/T) can- 
eigenvalue Aq. Equation ^ illustrates that even in the eels the thermal weight exp(-A„_,_p/T). By combining 
zero-temperature Umit it is in general not sufficient to Eq. dH]) as well as its five analogues with Eq. dHI), one 
consider only the corresponding ground state \z) = |0) can finally calculate the current by carrying out the re- 
when computing the thermal expectation value (. . .)o4i maining energy integrals numerically (or even analyti- 
Namely, it turns out that if |0) appears as any of the cally). More importantly, an intuitive understanding can 
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FIG. 2. (Color online) Gate voltage e dependence of the Josephson current (thick blue lines) as well as of the three lowest 
many-particle energies (thin red lines) for a single Anderson impurity coupled to BCS leads in the large-gap limit A = oo at 
phase difference (p/n = 0.5, different Coulomb interactions U and Zeeman fields B. The many-particle energies are shown in 
arbitrary units; those for finite B were shifted upwards for clarity as indicated by the arrow. The corresponding eigenstates 
are characterized by spin quantum numbers {s,Sz}- 



be obtained by interpreting J in terms of four successive 
tunnel events, each of which being described by one of 



the factors 
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III. RESULTS 
A. Single-level case 

In order to address the question in what respects the 
physics of our multi-level quantum dot in presence of 
a Zeeman field B Ki t can be interpreted in a simple 
single-impurity fashion, it is helpful to give a short re- 
minder of the latter. We focus on A = oo which allows 
to compute both the current as well as the many-particle 
eigenenergies/states. If one thinks of the physics being 
governed by some interplay of Kondo screening and the 
formation of Cooper pairs - and thus by the ratio IV/A 
- this large-gap limit might a priori seem highly non- 
generic. However, various prior studies^lii^ii^ii^ showed 
that treating A = cxd is sufficient to qualitatively un- 
derstand the single- level behavior (e.g., the parameter 
dependence of the phase boundary). In Section lilIC| we 
will illustrate that the same holds for our geometry. 

The effective Hamiltonian associated with the single- 
level case at large gaps can be obtained in complete anal- 
ogy with Eq. ([HDi^ii^i^ 
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= H 



dot ^2rcos((/./2) 
(e - C//2 4- B)n^ + (e 



C//2 - B)n^ + Un^n^ . 
(19) 



One can diagonalize it easily and eventually calculate the 
current from the phase derivative of the free energy [see 
Eq. ([Ti]) ]. The results are shown in Fig. [21 At zero Zee- 
man field i? = 0, we find that cither the ground state is 



a non-degenerate singlet with total spin {s = 0, s^ — 0} 
and the first excited state is a twofold degenerate dou- 
blet {s = 1/2, s^ = ±1/2} (the phase with J > 0), or 
vice versa (the tt phase where J<0; J = OatA=:oo). 
Decreasing U (or </>, or shifting e away from half filling) 
favors the phase; at a critical Uc (or 0c, or Ec), the sys- 
tem undergoes a first order level-crossing quantum phase 
transition as the nature of the ground states changes [see 
the B — Q curve J(e) in Fig.[2Ia); compare with Fig.[2Ib) 
where U < Uc\. The current is discontinuous and flips 
its sign. 

At finite B y^ 0, the degeneracy of the states with 
{1/2, ±1/2} is Ufted; but for small enough fields, the 
physics is still determined by an interplay of one non- 
degenerate and a pair of almost twofold degenerate 
states. As the singlet {0, 0} is not affected by the Zeeman 
field, the it phase becomes larger with increasing B; one 
can moreover trigger a 'Zeeman field O-tt transition' for 
U < Uchy applying a sufficiently large B [sec Fig. [H^b)]. 
The effect of B therefore seems similar to that of U; and 
even the lineshape of J(e, </>) at large B > and U = 
heuristically resembles the one at U > Uc and B = 0. 
However, the absence of an (almost) twofold degenerate 
state in the former case renders it meaningless to speak 
about an interplay of and tt behavior. 



B. Functional RG for generic parameters 

We start our discussion of the serial quantum dot by 
calculating the zero-temperature Josephson current J 
as a function of the gate voltage e for experimentally- 
motivated parameterai^ of the Coulomb interaction 
(C//r = 8), BCS gap (A/r = 2), and inter-dot cou- 
pling {t/r = 20). To this end, we solve the functional 
renormalization group flow equations ([5]) and ^ numer- 
ically and subsequently compute J from Eq. pi[) . The 
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FIG. 3. (Color online) Functional RG calculation for the gate-voltage e dependence of the Josephson current J flowing 
through the serial quantum dot for generic system parameters U /T — 8 (Coulomb interaction in units of the lead hybridization 
strength), A/F — 2, (jf>/7r = 0.5 (BCS gap and phase difference), i/F = 20 (inter-dot coupling), zero temperature, and 
various Zeeman fields: (a) B/T = 0.2, 2, 8, 16, 19.6 (from bottom to top at e = 0) (b) B/T = 19, 19.6, 19.8, 20.2, 20.4, 22 or 
(B - Bc)/F = -1, -0.4 - 0.2, 0.2, 0.4, 2 with Be = t; see the main text for details. 



results are shown in Fig. [H For small Zeeman fields, the 
lineshape of J(e) features two copies of the single-level 
current gate-voltage relation (compare, e.g., with Ref. [2J 
or with Fig. [2|) . This is easily understood from the fact 
that i?dot can be diagonalized by rotating to bonding and 
anti-bonding single-particle levels (see Sec. Ill D|) : they are 
well-separated at t S> F, and thus single-impurity behav- 
ior occurs if either is at the Fermi level (e « ±eo)- 

As B increases, the size of the regions with J < 
increases; but again it eventually becomes meaningless to 
strictly speak about or tt regimes. On the other hand, 
one readily notices that at B = Be = t)^ (i) The single- 
particle {U = 0) energies of the anti-bonding spin-up and 
bonding spin-down states become equalj^ they cross the 
Fermi level at e = 0. (ii) The smallest eigenvalue of 
Hdot (at U > 0) with two particles, which is the smallest 
overall one close to e = 0, becomes twofold degenerate; 
the associated spin configuration is either a singlet or 
(one out of a) triplet. Both motivates us to investigate 
whether for B « Be, e ~ the physics can again be 
described in a pure (and simple) single-impurity fashion. 
This turned out to be the case for normal leads where 
finite- -B Kondo ridges appear 4^ 

The Josephson current J(e) of the multi-level quan- 
tum dot with B Ki Be, e « indeed strikingly resembles 
the one for _B « 0, e « ±eo up to an irrelevant over- 
all sign [sec Fig. |3] and compare the curves in (a) which 
are connected by the arrow; compare (b) with Ref. |2J 
or with Fig. [5l^a)]. Most importantly, there are disconti- 
nuities associated with a sign change; and we find that 
the regimes of negative and positive current - from now 
on denoted and t:, respectively - feature precisely the 
same dependencies on system parameters as the and 
■K phases of a single impurity. Namely, decreasing U , e, 
\B — Bc\, or TT — (j> (an overall sign corresponds to a shift 
— >■ (f> + tt) favors the regime. Similar to the lineshape 
of J(e), the current-phase relation J{4>) around B sa B^ 



e w is analogous to the single-level case: It is half- 
sinusoidal (sinusoidal) in the (tt) regimei^ Those are 
all certainly rather heuristic arguments; but they can be 
fully supported by an analytic treatment of the large-gap 
limit (Sec. IIII Cp where one can access the exact many- 
particle eigenenergies/states. 

On the other hand, there is one clear and interesting 
difference to the single-level case. Namely, for U > Uc 
the current features another discontinuity at B — Be in 
addition to the O-tt transition taking place at U = Uc 
(or some e = £c, or (jj = c/jc)- It is associated with a 
change in magnitude of J but not a sign flip. This is illus- 
trated in Fig.[3{b); note that the curves for {B — Bc)/T = 
-0.4, -0.2 as well as for {B - Bc)/T = 0.2, 0.4 coincide, 
demonstrating that the current indeed exhibits a jump 
in the tt phase (a very detailed numerical investigation 
of this discontinuity seems superfluous as we can show 
analytically that it appears for A = oo; see Sec. IIIICp . 
On both sides - i.e. for all B ^ Be - of this suppos- 
edly first-order quantum phase transition, one observes 
TT phase behavior in the spirit outlined in the previous 
paragraph. Again, this is rather heuristic; and in order 
to obtain a more thorough understanding, we now analyt- 
ically investigate the atomic limit (A = oo) and moreover 
carry out a perturbation expansion in the tunnel coupling 
F. We finally note that this altogether scenario (summa- 
rized by Fig. [3|) is generic - it is not altered qualitatively 
if different local- and nearest neighbor interactions, level 
detunings or source-drain coupling asymmetries are in- 
troduced. 



C. Insights from the atomic limit 



As outlined in Section III CI our multi-level quantum 
dot geometry can be treated analytically in the large- 
gap limit even in presence of two-particle interactions; we 
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FIG. 4. (Color online) Exact results for the zero-temperature 
Josephson current in the infinite-gap limit (blue thick lines) 
compared to approximate FRG data obtained at A/F = 2000 
(orange dashed lines) for three different Zeeman fields. The 
parameters are the same as in Fig. [S] Thin red lines display 
the three lowest many-particle eigenenergies for B/F = 20.1 
(main panel) and B/F = 0.2 (inset) in arbitrary units. They 
can be characterized by total spin quantum numbers {5,5^}. 
For B Ri Be — t, there is a level-crossing phase transition 
between a (nearly) twofold degenerate ground state formed 
by {0, 0} and {1, —1} and a non-degenerate one as e is shifted 
away from the Fermi energy. This is the typical behavior 
of a single impurity at small Zeeman fields [which features 
an almost identical lineshape J(e); see the inset as well as 
Fig. [21^ a)]. In contrast to the latter, however, the magnitude of 
the current is different in each of the states {0, 0} and {1,-1}; 
there is another level-crossing phase transition at B = Be 
where the ground state changes from {0, 0} to {1, —1} or vice 
versa (this can be shown analytically; see also Fig. [5}. 



merely need to diagonalizc the effective Hamiltonian of 
Eq. (|12p in the underlying 16-dimensional Hilbert space 
and subsequently compute J from the phase derivative 
of the free energy [see Eq. (|14p ]. The current at A = 00 
is similar to the one at A/F = 2 (see Fig. |4] and com- 
pare with Fig. [3]) ; this is consistent with prior studies of 
the single-level casei^i^i^^ii^ In particular, its lineshape 
features single-impurity characteristics for B f^ Be = t, 
e sa accompanied by an additional discontinuity at 
B ~ Be in the tt regime. Again, J is larger a.t B < Be 
and even vanishes for _B > _Bc at A == 00. Before trying 
to understand this more thoroughly, let us just note that 
the results at A = 00 can be employed to benchmark our 
FRG approach. It turns out that despite the fact that the 
FRG embodies an approximation which can a priori be 
justified only for small values of the Coulomb interaction, 
its results compare decently with the exact reference even 
for fairly large U/T = 8 (compare thick solid and dashed 
lines in Fig. |4|) . This observation is again in agreement 
with the single-level case.— >^ 

Most importantly (and frequently referred to above), 
one can try to learn more about the physics by investi- 
gating the many-particle eigenvalues and eigenstates of 
ifeff and by characterizing them in terms of the square 
and z-component of the total spin s. For small B as well 



FIG. 5. (Color online) Atomic-limit (A — 00) Josephson cur- 
rent (thick line) and four lowest eigenenergies (thin lines; ar- 
bitrary units) as a function of the Zeeman field B at e = 
(other parameters are as in Fig.[3|. Since U > Uc, the ground 
state is almost twofold degenerate for B ^ Be, the involved 
states have a total spin {0,0} or {1, —1}, respectively. Thus, 
the spin configuration is a singlet for B < Be but (one out 
of a) triplet for B > Be - a,n additional level-crossing phase 
transition occurs at B = Be in the tt regime; the current is 
discontinuous. 



as e chosen such that either the bonding or anti-bonding 
level is close to the Fermi energy, one recovers (see the 
inset to Fig. 2]) the scenario outlined in Sec. IIII Al The 
ground state is either a non-degenerate spin singlet with 
{s = 0,Sz = 0} (the phase) or an almost twofold de- 
generate state associated with a free spin {1/2, ±1/2} 
(the IT phase; the ground state is exactly twofold degen- 
erate at B = 0). However, the very same degeneracy 
properties appear for B sa tt^ What we labeled and 
TT regime is always associated with a non-degenerate and 
a nearly (exactly a,t B = Be = t) doubly degenerate 
state, respectively. E.g., a.t U > Ue the ground state is 
nearly twofold degenerate around e sa but becomes non- 
degenerate as the gate voltage increases beyond ±ee (see 
the thin lines in the main panel of Fig. 3] and compare 
with the inset). At small U < Ue, the ground state is 
always non-degenerate whereas the first excited state is 
nearly twofold degenerate, and the current is continuous. 
Summarizing, this interplay of a non-degenerate and a 
(nearly) twofold degenerate state supports our conjecture 
of a level-crossing O-tt transition governing our multi- level 
quantum dot for Zeeman fields B ^ Be- 

The non-degenerate and pair of degenerate states in- 
volved arc characterized by total spin quantum num- 
bers {s — 1/2, Sz = —1/2} (0 phase) as well as {0,0}, 
{1,-1} (tt phase), respectively] '^ '^1^^ In the tt regime, the 
ground state is thus exactly twofold degenerate only at 
B = Be = t; for small deviations B sa Be, the degeneracy 
is lifted. The same holds in the single-level case; but there 
the current is obviously identical for both {1/2, ±1/2} 
(i.e., independent of the sign of B). For the problem 
at hand, however, the ground state spin configuration is 
different for B — Be < and i? — i3c > 0; it is a singlet 



0.04 






B/r=i8 


U/r=8, e=0 

t/r=20 


^^, 


- 


/B/r=2 


infinite A 


r^^^^^/'5\ 


PT 


X ^^XB/r=22 


— -FRG 


'^ B/r=30 -X,^ 




- ^^^ 


^.. 



10 



100 

A/r 



1000 



FIG. 6. (Color online) Zero-temperature Josephson current 
as a function of the BCS gap A for various Zeeman fields 
(all other parameters are as in Fig. [Sj. Solid lines were ob- 
tained by carrying out perturbation theory in the lead cou- 
pling r, dashed lines display FRG data. Dotted lines show 
the atomic-limit result. Remember that the current in the -k 
phase vanishes at A — >■ oo for B > Be but remains finite for 
B <Bc. 



in the former but (one out of a) a triplet in the latter 
case. This renders it reasonable that also J is different, 
and thus another level crossing phase transition occurs 
at B = Be in the tt regime (see Fig. [S]). Even though 
it might be intuitively clear that the current is smaller 
in the {1,-1} configuration^ - a triplet should prevent 
Cooper pair tunneling - we will now substantiate that by 
a perturbation theory approach. 



for B > Be). The e^-intcgrals in Eq. ((T6)) can then be 
carried out analytically: 



^A^oo 



F" sin 

2 V U+t+t 





U+t-t 



B<B, 

B>Br^ 



(20) 



For B < Be, one can moreover intuitively interpret every 
(of the few) non-vanishing terms in Eq. (|18p by asso- 
ciating a single tunnel event with each of the four fac- 
{n\dJl^\q). Fortunately, there are only two quali- 



tors 



tatively different scenarios: The two electrons occupying 
the do1(^ are first transferred to the right lead, and then 
the original dot configuration is restored by two electrons 
tunneling in from the left (see Fig. [T]); or alternatively, 
an additional Cooper pair from the left hops onto the dot 
and then further on to the right. In either case, there are 
two intermediate configurations that contain a single un- 
paired lead electron. Both give rise to a (quasiparticle) 
energy denominator E{ei) = \/i|"+"A^, but the current 
remains finite in the A — > c» limit due to the explicit 
prcfactor A^ in Eq. ((TB)) . On the other hand, it is quite 
clear that no similar process exists for B > Be where the 
dot is initially in a triplet state c^L c'L |vac) (i.e., occupied 
by two electrons with equal spin); but since there is no 
possibility for an intermediate state that does not exhibit 
a lead quasiparticle excitation, the current is suppressed 
more strongly at large gaps. This eventually provides a 
simple cartoon picture for the discontinuity at B = Be'. 
A Cooper pair can tunnel more easily through a singlet 
than a triplet dot configuration. 



D. Insights from perturbation theory 



IV. CONCLUSION 



The expansion of the Josephson current to fourth or- 
der in the level-lead coupling Hamiltonian Ht is given 
by Eq. (|16p . The imaginary time integrals can be carried 
out analytically [see Eq. (US])], but in general one needs to 
resort to numerics to perform the remaining two energy 
integrations. The resulting zero-temperature current as 
a function of the gap is shown in Fig. [HI One can easily 
understand that it again features a discontinuity around 
B = Be, e = 0; this is merely due to the fact that the 
smallest eigenvalue of Hdot corresponds to a spin singlet 
state (iLdL |vac) for B < Be but a triplet configuration 



'^li^bilvac 



' 1-'^'-^ for B > Ber^ and thus different terms con- 
tribute to Eq. ([T8| in each case. More generally, one ob- 
serves a decent quantitative agreement with FRG data 
even at A/F = 2 despite the fact that the expansion 
w.r.t. Ht can a priori be justified only for large gaps. 

Having established that perturbation theory gives sen- 
sible results, we now additionally turn to the limit of 
A — >■ cx); this large-gap case is again particularly simple 
since only a few contributions to Eq. ([T5|) as well as to 
its five analogues originating from different time order- 
ings remain finite for B < Be (and the current vanishes 



We have studied the equilibrium supercurrent J 
through a multi-level quantum dot in the molecular 
regime of large couplings t between the impurities. If 
a Zeeman field B ^ t is introduced such that the lowest 
two-particle energies (whose total spin configuration is 




FIG. 7. (Color online) The current through a singlet dot con- 
figuration interpreted in terms of four (1,2,3,4) successive tun- 
nel events in the basis of bonding and anti-bonding molecular 
levels. In the depicted contribution, only two intermediate 
states contain an unpaired lead electron (the one between 2 
and 3 does not). The latter require a quasiparticle excita- 
tion energy E{ei) = ^/ef + A^ and are thus unfavorable at 
large gaps. If the dot is initially in a triplet configuration 
dJj.dL |vac), no such process exists, and the current is sup- 
pressed more strongly. 
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FIG. 8. (Color online) Atomic-limit calculation of the Joseph- 
son current as a function of (a) the gate voltage e at ^/tt = 0.5 
and (b) the phase difference at e/F = 1.9 for U/T = 4, 
i/F = 20, B/F = 19.9 ('close to singlet-triplet degeneracy') 
and various temperatures T (J varies with increasing T as 
idicated by the arrows). Note that this is an exact result. 



either {s = 0,Sz = 0} or {s = l,Sz = —1}) of the iso- 
lated dot are equal, several characteristics of the system 
- such as the appearance of a O-tt transition as well as the 
corresponding lineshapes and parameter dependencies - 
are similar to the well-known case of a single impurity. 
On the other hand, the current features an additional 
discontinuity at B ~ Be = t in the n regime; it can 
be attributed to the fact that a Cooper pair can tunnel 
more easily through a spin singlet state than through a 
triplet one. This scenario should be generic for any more 
complex quantum dot geometry close a singlet-triplet de- 
generacy. 
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Appendix A: Finite temperatures 



The analytic treatment of the atomic limit outlined 
in Sec. Ill CI allows to readily compute J not only at 
T = but also for finite temperatures. This was not 
the main focus in the single-level case where the current 
at T = and A = oo always vanishes in the n regime 
[sec Fig. [21a)]; the discussion of the large-gap limit there- 
fore mostly concentrated on the parameter dependence of 
the phase boundary. For the problem at hand, J is non- 
zero in the tt phase for Zeeman fields B < B,. even at 
A = cx), and this motivates us to briefly present (exact!) 
results for the current at finite temperatures (see Fig. [8|) . 
Not unexpected, the discontinuities in J are smeared at 
T > 0; this is similar to the single- level case where non- 
zero r > were studied at finite gaps by virtue of the 
numerical renormalization group ^"'^"^ or quantum Monte 
Carlo32i2^ Moreover, the current-phase relation - which 
is generally (half-) sinusoidal in the tt (0) regimes - be- 
comes completely sinusoidal at sufficiently large temper- 
atures. This is once again completely analogous to the 
single-level case. 
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